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ABSOJRACrr 

Calculations  are  carried  out  to  determine  whether  images  of 
the  electronic  Fermi  surface  mi^t  he  observable  in  phonon  and  spin- 
wave  spectra  of  metaJ-s.  Semi -quantitative  estimates  are  made  of 
the  magnitude  of  image  effects  in  phonon  dispersion  curves  and  phonon 
lifetimes  in  typical  metals,  as  well  as  in  spin-wave  dispersion 
curves  in  rare-earth  metals.  These  estimates  lead  to  the  conclusion 
that  Fermi  surface  Images  may  be  experimentally  observable  in  favor¬ 
able  circumstances. 
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INTRODUCTION 

It  is  well  known  that  the  conduction  electrons  in  a  metal  have 
an  important  influence  on  the  effective  force  between  the  ion  cores. 
Similarly,  there  now  exists  good  evidence^  that  in  some  strongly  magnetic 
metals  the  effective  interaction  between  ionic  spins  is  in  large  measure 
due  to  the  mediation  of  the  conduction  electrons.  It  has  been  sug¬ 
gested  that,  as  a  result  of  these  roles  of  the  conduction  electrons, 
phonon  and  spin-wave  spectra  should  eadiibit  anomalies  which  reflect  with 
geometric  exactitude  the  shape  of  the  Fermi  surface.  These  emomalies 
are  referred  to  as  Images  of  the  Fermi  surface.  They  have  been  experi¬ 
mentally  observed  in  the  phonon  spectrum  of  lead  by  Brockhouse  et  ed.^ 

The  influence  of  the  conduction  electrons  on  the  effective  force 
between  two  ion  cores  can  be  described  in  a  sinple  way.  Displacement 
of  an  ion  from  its  equilibrium  position  virtually  excites  an  electjwn- 
hole  pair  through  coulomb  interaction,  and  the  pair  then  recombines-  by 
interaction  with  a  second  ion.  Calculation  of  the  energy  due  to  this 
process,  with  ion  positions  considered  as  parameters,  gives  an  effec¬ 
tive  potential  energy  of  interaction  between  two  ion  cores.  The  influence 
of  the  conduction  electrons  on  the  Interaction  between  ionic  spins  can  be 
found  similarly.  An  ion  is  able  to  excite  a  virtiml  electron-hole  pair 
by  coulomb  exchange  between  a  conduction  electron  and  one  of  the  ion's 
unpaired  electrons.  The  electron-hole  pair  recombines  by  interaction 
with  a  second  ion,  and  again  the  energy  due  to  this  process,  with  ionic 
spin  orientations  taken  as  pareiraeters,  gives  an  effective  Hamiltonian 
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for  interaction  between  ionic  spins. 

If  the  ions  are  displaced,  or  their  spins  are  oriented,  in  such  a 
way  th8,t  a  wave  of  wavenumber  £  propagates  throu^  the  lattice,  only 
electron-hole  pairs  of  wavenumber  c|;+K  will  be  created,  where  K  is  any 
vector  of  the  reciprocal  lattice.  As  the  energies  of  phonon  land  spin- 
wave  excitations  are  small  conpared  to  the  electronic  Fermi  energy, .  the 
energy  difference  between  the  states  of  the  whole  system  before  and 
after  virtual  creation  of  an  electrcn-hole  pair  is,  to  a  good  approxi¬ 
mation,  just  the  energy  of  the  electron-hole  pair.  For  small  enoxa^ 
q+K,  as  is  shown  in  Fig.  1  for  a  spherical  Fermi  surface,  it  is  possible 
for  electron-hole  pairs  of  zero  energy  to  be  created,  but  for  £+K  too 
large  the  lowest  energy  pair  that  can  be  created  has  finite  energy  pro¬ 
portional  (in  the  case  of  a  spherical  Fermi  surface)  to  |£tK|  ( |£+K|-2k5i). 

As  £  varies,  then,  zero  energy  pair  creation  abruptly  becomes 
impossible  on  surfaces  S  in  c^-space,  determined  by 
1)  3  +  K  =  2^ 

where  £  is  a  wave  vector  on  S  and  is  on  the  Fermi  surface.  The 
strength  of  the  effective  interaction  mediated  by  the  electrcaas  therefore 
sharply  decreases,  and  the  interaction  between  ions  changes  suddenly,  on 
the  surfaces  S.  For  phonons,  for  instance,  when  the  vector  ^(£+K)  passes 
with  increasing  £  from  inside  to  outside  the  Fermi  surface,  electron 
shielding  of  the  ions  will  be  sharply  reduced,  so  forces  between  the 
ions  will  become  stronger.  Then  the  phonon  frequency  will  sharply 
increase  on  the  surfaces  S  in  £-space.  On  the  other  hand,  if  i-(£+K) 
passes  from  outside  to  inside  the  Fermi  surface  with  increasing  £, 
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frequency  will  sharply  decrease.  Similar  anomalies  will  occur  ip  spin 
wave  frequencies. 

Because  real  pair  creation  cannot  occur  without  conservation  of 
energy,  it  is  apparent  that  the  lattice  abruptly  loses  the  ability  to 
give  energy  to  the  electron  gas  whenever  £  passes  throu^  one  of  the  svir- 
faces  S.  Thus  the  lifetimes  of  phonons  or  spin-waves,  considered  as 
functions  of  £,  should  show  sharp  breaks  on  the  surfaces  S. 

In  Part  I  we  calculate,  using  simple  models,  dispersion  curves  and 
lifetimes  for  phonons  in  a  typical  metal,  and  in  Part  II  we  calculate 
dispersion  curves  for  spin  waves  in  the  rare  earth  metals.  The  disper¬ 
sion  curve  anomalies,  are  shown  to  be  logarithmic  infinities  in  the 
derivative  of  frequency  with  £,  and  the  decay  rates  are  shown  to  display 
step-discontinuities.  The  magnitudes  of  these  effects  are  estimated. 

These  anomalies  may  in  favorable  circumstances  be-  observed  by  neutron 
or  X-ray  diffraction  measurements.  Their  observation  should  prove  very 
valuable,  as  the  location  of  the  anomalies  gives,  by  the  si]i5)le  construc¬ 
tion  of  Eq.  (l),  quite  direct  information  of  the  shape  of  the  Fermi 
surface . 

I.  IMAGES  OP  THE  FERMI  SURFACE  IN  PHONON  DISPERSION  CURVES  AND  PHONON 

LIFETIMES  IN  METALS 

Fermi  surface  images  in  phonon  dispersion  curves  for  lead  have  been  ' 
experimentally  observed,  as  mentioned  in  the  Introduction,  by  Brockhouse 
et  al.^  This  group  has  also  reported  that  no  image  effects  could  be 
detected  in  the  spectrum  of  sodium. 
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Here  we  attempt  to  make  semi-quantitative  estimates  of  the  size  of 
image  effects  in  phonon  spectra.  In  Sec.  A  phonon  dispersion  curves 
are  calciilated  for  a  metal  pictured  as  an  array  of  ions  immersed  in 
and  screened  by  a  uniform  gas  of  conduction  electrons.  The  magnitude 
of  the  calculated  image  effects  is  comparable  to  that  obseirved  in  lead.^ 
Sec.  B  refines  this  model  to  take  into  account  the  Bloch  character 
of  the  electron  wave  functions.  Using  assumptions  which  are  believed 
to'  be  valid  for  certain  matals,  including  sodium  and  aluminum,  we  find 

that  this  refinement  reduces  the  image  effects  by  a  large  factor,  of 

_o 

order  10  .  It  may  be  remarked  here  that  the  magnitude  of  image  effects 

is  directly  related  to  the  magnitude  of  certain  electron-phonon  Umklapp 
processes. 

In  Sec.  C  we  calculate  the  finite  phonon  lifetime  resulting  from 
decay  by  pair  excitation.  .  As  a  function  of  £  this  quantity  displays 
step  discontinuities  which  lie  on  the  same  surfaces  S  in  £  space  as  the 
frequency  anomalies.  These  discontinuities  may  be  experimentally  • 
observable  in  the  widths  of  phonon  frequencies. 

A.  The  Uniform  Electron  Gas  tfodel 

The  initial  step  in  our  calculation  of  dispersion  curves  is  the 
determination  of  an  effective  potential  of  Interaction  between  a  pair 
of  ions . 

We  first  consider,  following  the  calculation  of  Bardeen, ^  the 
electronic  charge  density  built  up  when  a  single  ion  is  placed  in  a 
uniform  electron  gas.  Let  the  ion  be  situated  at  the  coordinate  origin, 
and  denote  the  unscreened  ionic  potential  by  QV(r),  where  Q  is  the 
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ionic  charge.  Then  the  Hartree  Hamiltonian  for  an  individual  electron  is 

1.1)  H  =  Hg  +  i^V(r)  +  e^VgCr) 

where  Hq  is  the  kinetic  energy  of  the  electrons,  and  e^Vg(r)  is  the 
Hartree  potential  of  interaction  with  the  displaced  electronic  charg<!. 
This  potential  must  satisfy  Poisson's  equation 

1.2)  v2Vg(r)  =  -  4n  (N(r)  -  N^) 

where  N(r)  and  Nq  are,  respectively,  the  perturbed  and  unpertiurbed 
electron  densities.  We  introduce  the  voliime  of  the  crystal,  fi,  and 
Fourier  transfonns  according  to  the  notation 

1^3)  v(£)  =  J  V(r)  e'^E'^d^r. 

0 

Then,  by  first  order  perturbation  theory,  the  electron  wavefuncticns 
can  be  written  as 

1.4)  Y,  .  e 

where  the  Ej^  are  imperturbed  energies.  From  this  expression  n(£)  can  be 
c6ilculated.  Substitution  into  the  Fourier  transform  of  Eq.  (1,2/  gives 

1.5)  v^(e)  -  «(£)  . 

p2  +  8TTe2w(£) 


where 


1,6)  W(£)  =  -  n'  E 


k<kjp  E^  -  E^+^ 

and  the  sinmnation  nins  over  the  interior  of  the  Femri.  sphere  of  radius  kj, 


*T 


Because  of  the  charge  neutrality  of  the  entire  system,  the  zero 
wave  number  Fourier  conponents  of  all  potentials  are  consistently  dropped 
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This  expression  for  Vq(£)  can  be  substituted  into  Eq.  and  the 

electronic  charge  density  n(£)  confuted.  The  result  is 

1.7)  n(£)  =  2Qepgv(£)  W(jd) 

p^  +  8^reS^(£) 

6 

QMs  is  exactly  the  re sTilt .  found  by  Langer  and  Vosko. 

The  interaction  between  two  ions  can  now  be  viewed  as  having  two 
parts:  a  direct  interaction  between  the  two  ions,  and  an  interaction 
of  the  second  ion  with  the  electronic  charge  density  displaced  by  the 
first.  The  total  interaction  energy  when  the  second  ion  is  at  r  is  then 

1.8)  §(r)  = 

where 

1.9)  cp(e)  =  (^P^v^(£)Att)  I  1  -  8TTe^W(£) 

Y  p^  +  8TTe%(£) 

The  first  term  describes  direct  interaction  between  the  ions,  and  the 
second  the  reduction  due  to  electronic  screening. 

Calculation  of  the  vibration  frequencies  of  a  Bravais  lattice  is  now 
sinple.  We  denote  the  equilibriinn  position  of  the  m'th  ion  by  and 
its  excursion  by  Then,  to  first  order  in  the  excursions,  the  force 
on  the  m'th  ion  is 

1.10)  i£iji  =  m  ^mn  ~ 

=  e^-’^-“'-“^£-{£in-£n) 

1 

where  goes  over  the  N  ions  of  the  lattice,  excluding  m=n.  IBie  equations 
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of  motion  for  the  ions  are 

1,11)  =  M  r„ 

•’  '  — m  — m 

where  M  is  the  ionic  mass.  Eqs.  (l>10)  and  (l,ll)  are  solved  by  the 
AXXB&tZ 


1.12)  rj„  =  e  cos(_a*Rjn  -  cot) 

where  e  is  a  unit  polarization  vector,  £  the  propagation  vector,  and 

1.13)  u)2  =  (n/M  fi)  j^A(^)  +  S^q(a(K+^)  -  A(K))j 


The  sum  runs  over  all  vectors  of  the  reciprocal  lattice,  and 
1,14)  A(£)  =  (^-e)^  cp(p)  = 


=  (^PAtt)  (£-e)^  P^v^(£) 


8TTe^W(£) 
p^  +  8TTe^W(£) 


The  Fermi  surface  images  in  the  phonon  dispersion  curves  arise 
from  the  character  of  the  function  w(£).  A  straightforward  integration 
of  Eq.  (l,6)  yields  for  W(£) 


1.15)  W(£)  =  (m%/8TT2fi2)  In 

A  plot  of  W(£)  is  presented  in  Fig.  2.  It  should  be  noted  that  the 
gradient  of  W(£)  has  a  logarithmic  infinity  at  p=2kp.  ihe  phonon  fre¬ 
quency  as  defined  by  Eqs.  (l,13)  and  (l,l4)  then  has  logarithmically 
infinite  gradient  with  respect  to  £  on  the  surfaces  S  in  £  space  defined 
by 

1.16)  \K+g\  =  2kp 

where  K  is  any  vector  of  the  reciprocal  lattice.  Consequently,  "kinks" 
occur  in  the  dispersion  curves  on  these  surfaces,  (in  the  more  general 
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case  of  a  non-spherical  Fermi  surface,  Eq.  (I,l6)  shoixld  "be  repljiced  by 
Eq.  (l)  of  the  introduction.  Hie  case  of  a  non-spherical  Fermi  surface 
is  treated  somewhat  more  fully  in  Appendix  1. ) 

The  dispersion  cuirves  (l,13)  have  been  computed^  for  constants 
corresponding  to  aluminum,  for  propagation  vectors  in  directions  of  high 
symmetry.  Dispersion  curves  for  £  in  the  (ill),  direction  of  reciproceJ. 
space,  for  both  directions  of  polarization,  are  plotted  in  Fig.  3* 

While  the  model  of  this  section  almost  certainly  overestimates  the 
prominence  of  image  effects  in  aiuminum,  (see  Sec.  B)  the  qualitative 
similarity  of  Fig.  3  to  the  observed  phonon  spectrum  in  lead  (Fig.  4)  is 
quite  striking. 

B.  Bloch  Electrons  in  the  Bardeen  Approximation 

In  a  classic  paper  on  electrical  conductivity  of  metals,  Bardeen^ 
succeeded  in  treating  the  interaction  of  lattice  vibrations  with  Bloch 
electrons  by  making  certain  simplifying  assunptions  about  the  electron 
wavefunctions .  In  this  section  we  estimate  the  magnitude  of  image 
effects  under  the  same  assumptions. 

We  first  calculate  the  change  in  electronic  charge  density  set  t;® 
when  one  ion  of  the  lattice  is  displaced  from  its  equilibrium  position. 

We  take  as  electron  Hamiltonian 
1,17)  H  =  Ho  +  ^  rj^.  VV(r)  +  e^Ve(r) 

where  Hq  contains  the  electron  kinetic  energy,  the  energy  of  interaction 
with  the  ions  in  their  equilibrium  positions,  and  the  Hartree  Interaction 
with  the  other  Bloch  electrons;  displacement  of  ion  1  from  its 
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equilibrium  position  at  the  origin;  ^V(r)  is  the  ionic  potential  and 
e^Vg(r)  is  the  self-consistent  electron  potential  due  to  displaceitent 
of  the  ion. 

We  first  assume,  following  Bardeen, 5  that  the  matrix  element 
(k+£|vv(r) |k),  taken  between  Bloch  states,  has  negligible  k  dependence. 

(This  is  certainly  so  if  the  Bloch  functions  can  be  written  in  the  form 
UQ(r)  ei^’£  .)  We  then  define  v*(£)  according  to  the  equation 

1.18)  £  v*(£)  s  -i  n  (k+£|vv(r)  |k) 

At  this  point,  to  avoid  making  a  fully  self-consistent  calculation 
with  the  Bloch  functions,  we  recognize  that  the  important  effect  of  the 
term  e^Vg(r)  is  the  screening  of  the  ionic  potential,  and  we  follow  Betr- 
deen^  in  saying  that  screening  by  HLoch  electrons  will  be  roughly  the 
same  as  free-electron  screening.  Then,  instead  of  keeping  the  term 
e^Vg(r)  in  Eq.  (l,17),  we  simply  reduce  v*(£)  by  the  factor  p^/(p^+8ne%(£)). 
Since  the  effect  of  the  screening  is  large  only  for  small  £  while  our 
region  of  interest  is  p^k^,  we  do  not  expect  this  approximation  to 
affect  our  results  seriously. 

With  these  assumptions  the  electron  wavefuncticns,  perturbed  by  the 
displacement  of  ion  1,  become 

1.19)  \  y*(t)  <P 

where  c(^  are  the  Bloch,  functions  and  Efc  are  the  eigenvalues  of  Hq. 

A  second  displaced  ion  interacts  electrostatically  with  the  first 
ion  and  the  displaced  electron  charge.  The  energy  of  intera,ction  with 
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the  second  ion  at  Rg+£2 


:,20)  ^  i:  e‘-E-.tel-S2)(e‘-E-<il-£2>  p272(E)-(E-ri)(E-l2)^^!^S>^^ 


+  8rre%(£) 


where  again  the  first  term  describes  direct  interaction  between  thb  ions^ 
and  the  second  term  the  reduction  due  to  electron  screening. 

The  calculation  of  vibration  frequencies  proceeds  essentlELlly  e^s  ih 
Sec.  A,  and  we  find 

1,21)  u)^  =  (n/M  d)  jA(q)  +  E  (  A{K+3)  -  A(K) 

where 


1,22)  A(e)  -  ^  (p.e)2  r pV(e)  -  P^v*^(e)  (i  .  5  ) 

p^  +  8TTe^W(£) 

Except  for  substitution  of  the  factor  v*^(£)(l  -  6^  for  v^(£),  this 


res\ilt  is  identical  to  the  result  for  the  uniform  electron  gas  model. 
The  size  of  image  effects  is  determined  by  the  magnitude  of  v  (p)  for 
p=2kji.  For  sodium,  Bardeen  gives ^ 

1,23)  ^v*(p)  =  (C  +4TT^/p2)  G(prs) 


where  rg  is  the  radius  of  the  unit  cell  approximated  as  a  sphere,  and 
1,24)  G(x)  =  3x  sin  X  -  x  cos  x  )  G(0)  =  1 

The  constant  C  is  related  to  the  curvature  of  the  Bloch  functions  near 
the  edge  of  the  vinit  cell,  and  is  small  (less  than  ten  percent  of  4TTe2/p2 

Q  . 

in  sodium  for  p=2kj,. )  Here  we  ignore  C  and  note  that  the  in^jortant 
modification  of  the  result  of  the  last  section  is  a  reduction  of  electron 
shielding  terms  by  the  factor  G^(pr  ).  Now  G^(x)  is  a  rapidly  decreasing 

■S 
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function  of  x.  In  fact,  for  all  metals  with  ^  in  the  first 

Brillouin  zone  and  k/o,  is  small  con^jared  with  unity,  in  general  of 

order  10“^.  As  a  resTiLt,  the  electronic  contribution  to  A(£)  (the 

second  term  in  Eg.  (l,22))  can  be  to  a  good  approximation  neglected  in 

the  occurring  in  Eg.  (l,2l).  In  particular,  the  magnitude  of  the 

-P 

image  effects  is  reduced  by  a  factor  of  order  10  . 

These  considerations  are  in  the  same  spirit  as  the  discussion  of 

O 

the  elastic  constants  of  sodium  by  Bardeen  and  Pines,  who  obtain  agree¬ 
ment  with  e3g)eriment  to  within  ten  percent.  Their  results  can  be 
obtained  from  Egs.  (l,2l)  and  (l,22).  if  the  second  term  in  (l,22)  is 
dropped  for 

We  expect  this  estimate  of  Bloch  wave  effects  to  be  someidiat  drastic. 
It  is  probably  realistic  for  sodium  and  for  certain  other  metals  (includ¬ 
ing  al\uiiinum)  whose  characteristics  are  more  fully  discussed  in  J^apendlx 
2.  It  should  be  noted  that  image  effects  are  small  in  these  cases  for 
the  reason  that  Umklapp  processes  do  not  mate  an  tm^rtant  contribution 
to  the  electron-ion  interaction.  This  should  not  be  true  in  general. 

If,  for  some  metals,  Umklapp  processes  do  make  an  iaportant  contribution, 
we  shoxild  ejcpect  this  fact  to  be  reflected  in  some  characteristics  of 
the  metal.  Thus,  clues  may  appear,  from  an  examination  of  metallic 
characteristics,  as  to  whether  Fermi  surface  images  are  observable,  no 
such  clues  sxiggest  themselves.  First,  resistivity  of  a  metal  with 
ini)ortant  Umklapp  processes  shoiald  be  someidiat  hig^r  than  is  to  be 
expected  from  simple  scaling.  Second,  since  only  Uniklapp  processes 
screen,  the  treinsverse  lattice  vibrations,  metals  in  idil<di  Ihikl^pp 
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processes  are  negligible  will  have  transverse  frequencies  entirely 
determined  by  direct  ionic  interaction.  We  have,  of  course,  neglected 
ion-overlap  contributions  to  the  frequency;  so  that  transverse  frequencies 
will  be  higher,  by  some  unknown  amoiuat,  than  our  foinnulas  would  predict; 
but  at  least  the  existence  of  unusually  low  transverse  frequencies  would 
imply  that  Uipklapp  processes  were  contributing.^  In  both  of  these  char¬ 
acteristics,  lead  shows  the  properties  of  a  metal  with  substantial 
Umklapp  effects:  the  resistance  is  unusually  hi^  and  the  transverse 
frequencies,  as  judged  by  the  Debye  temperature,  are  unusually  low.  On 
the  other  hand,  aluminum  and  sodium  have  resistance  and  Debye  tenpera^ 
ture  in  line  with  the  usual  run  of  metals. 

In  summary,  the  non-observation  of  image  effects  in  sodliun  and 
aliunlnTam  is  consistent  with  other  enpirical  properties  of  these  metals,  ■ 
and  can  also  be  concluded  from  a  priori  calculations.'  On  the  other  hand, 
the  obser-vation  of  such  effects  in  lead,  while  consistent  with  other 
empirical  information,  has  not  yet  been  explained  on  the  basis  of  a 
de'tailed  calculation. 

C.  Phonon  Lifetime 

Several  decay  mechanisms  can  contribute  to  shortening  of  the  phonon 
lifetime,  including  interactions  with  other  phonons  and  various  processes 
which  transfer  energy  to  the  electron  gas.  We  here  discuss  the  decay 
process  which  gives  sinple  information  about  the  shape  of  the  Fermi 
surface;  that  is,  decay  of  a  phonon  into  a  single  electron-hole  pair.. 

In  Appendix  3  decay  ra't^s  for  inpor'tant  conpeting  decay  processes  are 
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estimated,  and  it  is  found  that  the  single  pair  creation  discussed  here 
gives  the  most  important  contribution  to  phonon  decay. 

We  find  the  lifetime  of  a  phonon  exposed  to  the  hazard  of  decay  liito 
an  electron  hole  pair  by  calculating  the  rate  of  energy  absorption  by  the 
electron  gas  when  the  lattice  is  driven  in  one  of  its  normal  modes.  We 
assume  that  the  electrons  see  screened  ions,  so,  the  electron- ion  matrix 
elenfent  is  taken  to  be  reduced,  as  in  Sec.  B,  by  the  factor  p^/(p^+8TTe%(£)). 
From  first  order  perturbation  theory,  we  find  for  the  probability  per  unit 
time  that  the  state  |k+£)  be  excited  from  the  state  |k),  by  Interaction 
with  an  array  of  ions  vibrating  as  in  ( I, 12) 

^  2  2/  2  2  (£-e)^ 

I>25)  ®  - 


^£,K+£  ^ 


p'^  +  8TTe^W(£) 

X  («(%+£  -  Ek  +  *U))) 

The  total  energy  absorbed  per  unit  time  by  the  electron  gas  is 

1.26)  E'  =  W^+£,k  (®k+£  "  %)  ■ 

where  runs  over  all  states  below  the  Fermi  surface,  and  Z  runs 

K<;Kp  £ 

over  all  electron  states.  is,  except  for  the  energy  6-functions, 

independent  of  k.  Converting  iti'fco  an  integration  and  doing  the 

integration  over  the  6 -functions  gives  a  factor 

1.27) 


m*2(D2/ 2TT2ft2p 
0 


for  p<2kj. 
for  p>2kp 


!nie  total  energy  absorption  per  unit  time,  expressed  as  a  fraction  of 
the  total  vibrational  energy  of  the  lattice,  can  be  Identified  as  the 


phonon  decay  rate,  or  Inverse  lifetime* 
1,28)  1/t  =  E  B(K+^) 

where 


1.29)  B(£)  = 


(£.e)2  p3v*^(p)/(p2  +  8ne%(p)) 


for  jxSkj. 


0  for  p>2kp 

The  decay  rates  confuted  from  Eq.s.  (l,28)  and  (l,29)  display  step- 
discontinuities  on  surfaces  in  £  space  defined  by  |£+K|  =  2kp.  These 
decay  rates  have  been  coniputed,  taking  for  v*(p)  the  nearly- free- 
electron  matrix  element,  for  constants  corresponding  to  aluminum.  Curves 
for  both  directions  of  polarization  for  £  in  the  (ill)  direction  of 
reciprocal  space  are  plotted  in  Pig.  5* 

Again  it  should  be  noted  that,  if  v*(£)  falls  off  faster  than  the 
nearly-free-electron  approximation  predicts,  the  magnitude  of  the  dis- 
continuities  will  be  correspondingly  reduced.  In  particular,  for  the 
case  of  sodium  and  aluminum,  only  the  K=0  term  should  appreciably  con¬ 
tribute  in  Eq.  (l,28),  and  no  significant  breaks  should  appear. 


II.  IMACaiS  OF  THE  FERMI  SURFACE  IN  SPIN-WAVE  DISPERSION  CURVES  AND  SPIN- 
WAVE  LIFETIMES  IN  RARE  EARTH  METALS 

The  rare  earth  metals  from  gadolinium  to  thulium  are  expected  to 
be  good  subjects  for  detection  of  Fermi  svirface  images  in  spin-wave 


*Neglecting  Umklapp  processes,  this  result  agrees  with  Pippard's^^ 
result  for  attenuation  of  ultrasonic  waves,  in  the  limit  where  £  is  small, 
thou^  not  so  small  as  the  reciprocal  electron  mean  free  path  from 
impurity  scattering. 
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dispersion  curves.  This  expectation  arises  "because,  in  view  of  the 
smallness  of  the  4f  shell  in  these  metals,  overlap  "between  4f  shells  of 
nei^boring  ions  should  be  quite  small.  Mie  ionic  spins  can  then  be 
expected  to  be  coupled  primarily  throu^  exchange  scattering  of  con¬ 
duction  electrons. 

In  this  section  we  develop  an  effective  Hamiltonian  for  interaction 
between  ionic  spins,  first  choosing  suitable  ionic  potentials  for  inter¬ 
action  between  conduction  electrons  and  ionic  spins.  The  change  in 
total  energy  of  the  system  due  to  these  interactions  is  then  found  by  a 
perturbation  calculation  carried  out  to  second  order,  paralleling  the 
ceilculation  of  de  Gennes^^  for  Curie  temperatures  of  the  rare  earths 
and  the  calculation  of  Ruderman  and  Kittel^  for  the  coupling  of  nuclear 
spins.  The  effective  Hamiltonian  for  interaction  of  ionic  spins  takes 

the  form  of  a  sum  of  Heisenberg  scalar  product  interactions.  The  spin- 

14 

wave  excitations  are  foixnd  by  usual  methods. 

Approximate  numerical  calculations  of  the  spin-wave  frequencies  as 
functions  of  wave  vector  are  carried  out  for  directions  of  high  symmetry 

I 

so  that  magnitudes  of  the  anomalies  can  be  estimated.  Uliroughout  the  dis¬ 
cussion  it  is  assumed  that  the  ground  state  of  the  ion  spin  system  is 
ferromagnetic.  The  extension  of  our  results  to  the  antiferromagnetic 
case  is  taken  up,  briefly,  in  Appendix  4. 

Because  of  the  small  radius  of  the  4f  shell,  we  take  the  exchange 
interaction  between  conduction  electrons  and  ionic  spin  to  have  6 -function 
spatial  dependence.  The  single  electron  Hamiltonian  can  then  be  written 
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II,l)  H  =  GE  6(R„  -  r)  J^-S 

’  '  m  m  —  — m  —  “  , 

where  Jjn  and  S  are,  respectively,  the  angular  momentinn  operator  for  the 
m'th  Ion  and  the  electron  spin  operator.  is  the  displacement  of  the 

electron  relative  to  the  position of  the  m'th  ionj  the  sum  runs  over 
H.1 1  N  ions  of  the  lattice.  The  scalar  product  form  Jjq-S  has  been  discussed 
by  Kasuya^^  and  by  Liu. rphg  interaction  (ll,l)  has  been  postixlated  by 
de' Gennes  and  Friedel^  in  their  ce^culation  of  spin-disorder  resistivity 
in  gadolinium;  they  were  able  to  deduce  an  approximate  value  for  the 
coupling  constant  G  from  experimental  resistivities.  De  Gennes^^  has 
suggested  that,  for  the  metaJ-s  from  gadolinium  to  thulium,  G  can  be  taken 
proportional  to  (g-l)j,  where  g  is  the  Lande  factor. 

The  effective  Hamiltonian  for  the  ionic  spins  is  fOTind  by  calcu¬ 
lating  the  energy  of  the  electron  gas  in  a  given  configuration  of  ionic 
spins.  The  part  of  this  energy  which  couples  the  spins  of  different 
ions  is,  to  lowest  order  in  the  coupling  constant  G, 

n,2)  X 

X  (k,s|E^6(Rj^-r)  J^-S  |k+£,s’)  (k+E,s'|E^  6(5^-E)  i„-  S|k,B) 

^ere  over  all  electron  states,  characterized  by  momentum 

eind  z-component  of  spin  s,  which  lie  beneath  the  Fermi  surface,  and 

runs  over  all  electron  states.  S'  runs  over  the  N  ions  of  the  lattice, 

n 

excluding  m=n.  The  energies  Ejj  have  been  taken  to  be  spin  Independent. 
Performing  the  spin  sums  and  evaluating  the  matrix  elements  gives 
11,3)  Heff  =  2m  4  Jm'Jn 


s 


:  .  18 
j  '  F(Rjjj-Rjj)  is  the  long  range  oscillating  intera,ctioh  introduced  "by  Ruder- 

i  man  and  Kittel,^  and.  is  given  by 

j  II, W  t(r)  . 

j’  where  n  is  the  volume  of  the  crystal. 

j  We  now  substitute  into  Eq.  (ll, 3)  the  Hols tein-Primakoff  operators 

;  defined  by^^ 

!  11,5)  (2J)i(l  -  aJa^SJ)^  = ‘^mx  +  ^^11^ 

.  (2J)i4(l  ■  4%/2J)^  =  Jmx  -  iJnor 

J  -  a^in  =  Jmz 

GSie  resulting  expression  can  be  expanded  to  give  terms  of  O'th,  2nd,  4ti]i, 

>  ...  order  in  the  field  operators.  We  omit  terms  of  hi^er  than  second 

I.  order.  (The  higher  order  terms  represent  interactions  between  spin 

^  ' 

waves,  and  give  no  contribution  when  only  one  spin  wave  is  present  in 
'  the  lattice.)  Then 

11,6)  F(Rj^-Rq)  j^l  +  {2/j)(-a^  +  a^jj 

The  hexagonal-close-packed  structure  of  the  rare  earths  requires 
that  a  distinction  be  made  between  the  two  types  of  lattice  site.  ,  We 
regard  the  ions  as  lying  on  two  hexagonal  Bravais  lattices,  whose  rela¬ 
tive  displacement  is  given  by  the  vector  D.  The  position  of  an  ion  on 
the  first  lattice  will  continue  to  be  written  as  m=  1,  2,  . . . ,  -gNj 
and  the  corresponding  spin  operators  will  be  written  as 
For  the  second  lattice  the  position  and  spin  operators  will  be  Rjjj+D  and 
We  substitute  into  Eq.  (ll,6)  the  transformations 
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II, T)  =■  (2/H)i 

a,'®’  .  (2/N)i  el3-(S„+2>,^(°) 

Here  the  sums  run  over  all  wave  vectors  £  in  the  first  Brillouih  zone  of 
the  hexagonal  Bravais  lattice.  Omitting  the  constant  term  in  Eq.  (ll^6) 
we  get 

11.8)  Hgff  =  g2je^  ^ 

where 

11.9)  \  (f(R)  +  F  (R+D)  -  P(R)  e^^'S) 

B^  =  F(R+d)  e^9.*  (£+£) 

Here  the  sums  run  over  all  the  displacements  R  of  a  single  BraVals  lattice, 
including  R^=0. 

Finally,  to  diagonalize  the  Hamiltonian  given  hy  Eq.  (II,8),  we  Ujse 


11,10)  h^  =  2'2 

1 

II 

The  re  stilt  is 

11.11)  Hgff  =  g2jS^ 

For  simplicity,  our  computations  will  he  limited  to  the  "acoustical"  spin 
waves,  whose  fre^encies  are  given  hy 

11.12)  iiu)  =  G%  (A^  -  .|b^|). 

To  carry  out  the  computations,  we  need  the  Fourier  transform  of  F(R), 
which,  hy  inspection  of  the  definition  (ll,4)  of  F(r),  is  just 

11.13)  fl  ^  ri5;<kp  =  ^(£) 


a*(o)a  ( 


O)). 


£  £  £ 


£  £  £ 
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where  W(£)  is  just  the  function  encountered  in  the  treatment  of  phonon 
dispersion  curves,  (plotted  in  Fig.  2)  which  has  logarithmically  infinite 
gradient  at  p=2kp>.  When  the  lattice  sums  indicated  in  Eq..  (ll,9)  are 
carried  out,  we  find 

II, li^)  =  i(N/n)  (w(K+£)  -  W(K)  -  w(K)cos(K.D)^ 

=  i(N/n)  %  W(K+2)  e^-*^ 


with  the  sums  running  over  all  vectors  K  of  the  reciprocal  lattice. 

The  spin  wave  frequencies  are  obtained  by  substituting  Eq.  (ll,l4) 
into  Eq.  (ll,12).  The  reciprocal  vectors  form  a  hexagonal  lattice.  For 
wave  vectors  pointing  in  the  direction  of  the  c-axis  in  reciprocal  space, 
the  formula  for  the  acoustical  frequencies  takes  the  particularly  simple 
form 


11,15)  =  G%(E/2n)%  ^W(K+ci)  -  W(K)^  (1  -  cos(K.D)^ 

Approximate  numerical  computations  of  the  dispersion  curve  in  the 
directions  of  the  c  and  a  axes  have  been  carried  out.  The  results,  are 
.  shown  in  Fig.  6.  Discrete  sums  over  the  first  3OT  points  in  reciprocal 
space  were  taken.  The  siun  ^W(K+q)  -  W(K)  ^  e^— *£  was  foimd  to  converge 
sufficiently  rapidly  that  more  distant  points  co’ild  be  neglected.  '  The 
svan  Ejj  ^w(K+^)  -  W(K)^  over  the  more  distant  points  was  approximated  by 
a  suitable  integral. 

As  in  the  phonon  ease,  the  frequencies  defined  by  Eqs.  (ll,12)  and 
(ll,li)-)  have  logarithmically  infinite  gradient  with  respect  to  £  on  the 
siirfaces  S  in  £  space  defined  by  Eq.  (I,l6)  or  Eq.  (l)  of  the  Introduction. 
Consequently,  "kinks"  occur  in  the  dispersion  curves  on  these  surfaces. 
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From  the  computed  curve,  these  kinks  can  be  seen  to  be  of  the  order  of 
two  percent  of  the  maximum  spin  wave  frequencies.  They  appear,  there¬ 
fore,  to  be  within  the  range  of  experimental  detection. 

It  will  be  observed  that  the  frequencies  become  negative  for  £  in 
the  c  direction.  The  caJ-ciiLations  have  the  feature  that  the  frequencies 
are  extremely  sensitive,  both  in  magnitude  and  sign,  to  small  changes  in 
the  shape  of  the  Fermi  surface.  However,  the  magnitude  of  the  kinks  is 
not  seriously  affected  by  such  changes.  The  implications  of  negative 
dispersion  curves  are  discussed  more  fully  in  Appendix  ll-. 

A  few  remarks  about  the  effects  of  finite  ion  size  seem  in  order. 
This  effect  is  extremely  difficult  to  treat  quantitatively.  However,  it 
should  be  noticed  that  to  some  extent  the  effect  is  allowed  for  in  the 
determination  of  G  from  observed  resistivities,  and  that  the  dimension¬ 
less  parameter  that  measures  the  iinportance  of  the  effect,  namely  kpr^, 
where  Tq  is  the  "radius"  of  the  4f  shell,  has  the  approximate  value  0.4. 
It  is,  therefore,  not  to  be  expected  that  our  estimate  of  the  relative 
size  of  Fermi  surface  images  is  seriously  in  error. 

It  should  be  noted  that  the  appearance  of  FermL  sTurface  images  in 
the  spin-wave  dispersion  curves  is  in  no  way  dependent  on  the  special 
features  of  the  rare  earth  metals.  In  general,  the  conduction  electrons 
of  a  metal  can  be  expected  to  take  some  part  in  the  interaction  between 
ionic  spins,  so  "kinks"  should  appear  in  the  spin-wave  dispersion  curves. 
Their  relative  size  will,  however,  be  reduced  as  the  electrons  play  a 
less  important  part  in  the  total  interaction  between  ionic  spins. 
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Fermi  surface  image  effects  in  the  spin-wave  lifetime  will  be  men¬ 
tioned  only  briefly  here.  The  contribution  to  spin-wave  decay  rate  of 
the  process  (spin-wave  -*  electron-hole  pair)  can  be  evalmted  strai^t- 
forwardly,  in  direct  analogy  to  the  phonon  case,  and  will  give  a  decay 
rate  which  has  step-discontinuities  on  the  surfaces  S  in  _^-space.  Other 
contributions  to  the  spin-wave  decay  rate  are,  however,  not  well  under¬ 
stood,^^  so  the  relative  size  of  the  FerniL  surface  images  is  difficult 


to  estimate . 
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APPENDIX  1 

Connection  Between  Anomalies  of  Frequency  and  Decay  Rate 

We  investigate  in  this  appendix  the  relationship  between  phonon 
frequency  and  phonon  decay  rate,  in  order  to  clarify  the  connection 
between  the  two  sorts  of  spectral  anomaly.  We  shall  find  that  a  close 
relationship  exists,  connecting  the  nature  and  magnitude  bf  Fermi  surface 
images  in  the  phonon  dispersion  curves  with  the  nature  and  magnitude  of 
Fermi  surface  images  in  the  phonon  decay  rate  curves. 

In  the  dispersion  curve  calculations  pf  Part  I,  it  was  assumed 
that  ion  motions  were  slow  enou^  that  screening  electrons  were  able  to 
follow  a  moving  ion  perfectly.  In  Sec.  a  of  this  appendix  we  relax  this 
assinnption,  calculating  the  electron  charge  distribution  time-dependently ' 
in  a  lattice  of  moving  ions.  The  force  exerted  on  a  given  ion  by  this 
electron  charge  distribution  is  found  to  be  not  altogether  in  phase 
with  the  position  of  the  ion,  so  the  ion  does  net  work  on  the  electron, 
gas,  and  there  is  net  transfer  of  energy  into  the  electron  gas  from  the 
vibrating  lattice.  The  rate  of  energy  absorption  can  be  related  to  the 
phonon  lifetime;  the  phonon  decay  rate  curves  and  dispersion  curves  are 
then  found  to  be  given  by  real  arid  imaginary  parts  of  the  same  function. 

In  Sec.  b  the  properties  of  the  function  which  gives  dispersion 
curves  and  decay  rates  are  investigated  for  coii5)lex  frequency  ou.  It  is 
shown  that  "kinks"  in  the  dispersion  curve  and  breaks  in  the  decay  rate 
arise  from  the  behavior  of  this  function  near  a  branch  cut  on  the  real 
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axis  of  the  cor^jlex  cu-plane.  The  relative  magnitudes  of  the  anomalies 
are  con^iared. 

In  Sec.  c  the  behavior  of  the  phonon  dispersion  curves  and  phonon' 
lifetimes  is  considered  for  the  case  of  a  non-spherical  Fermi  svirface. 

It  is  shown  that  Images  of  convex  parts  of  the  Fermi  surface  take  the 
same  forms  as  the  images  which  occur  in  the  case  of  a  spherical  Fermi 
surface.  The  relative  magnitudes  of  freqoaency  and  decay  rate  anomalies 
afe,  however,  found  to  be  somewhat  changed. 

a.  Calculation  of  Frequency  and  Decay  Rate 

We  proceed  to  calculate  the  time -dependent  electron  charge  distri¬ 
bution  in  a  lattice  ;diere  the  ionic  excursions  are  given. *by 

Al)  r  =  e  cos(q.R  -ojt) 

— m  *-  -*  — m 

where  Rjjj  is  the  equilibrium  position  of  the  m'th  ion.  Bie  electrons  are 
assumed  not  to  Interact  with  each  other,  though  the  largest  effect  of 
electron-electron  interaction  is  taken  into  account  by  assuming  that  the 
electrons  see  screened  ions.  An  electron  is  thus  assumed  to  see  a  poten¬ 
tial  (e^r)  from  an  ion  located  at  the  origin.  The  constant  a  is 

taken  to  be  the  Thomas -Fermi  screening  parameter  (2kp/TTaQ)5^,  where  a^ 
is  the  Bohr  radius.^ 

The  time  dependent .  electron  wave  functions  can  be  written 

A2)  fj^(r,t)  =  (e^Ck.r-E^t/fi]  ^ 

where,  from  perturbation  theory 

-^3)  s-k+£,k^^)  =  -(i/ft)  J  e^^’e^(%-t2"®k)'<^/*V:k-*£l^^V(r-RjQ-rm)|l^dt' 


"W: 


r 
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lEhe  screened  ion  potentials  V(r)  have  "been  taken  to  he  turned  on  slowly 
over  a  long  period  Of  time  hy  the  factor  e  '  .  T)  will  he  taken  small, 

so  that  Tlt«l  for  all  times  t  of  interest.  The  matrix  element  gives 
Alf)  <k+E|Qe£^V(r-R„-rJ|  k>  = 

r  I  ,  >  /  (»+'/'  .  ,i(«t 


+  „  -e- 

£,-K-E 


')] 


and  after  doing  the  time  integrations  the  electron  charge  density  can 
he  found  to  he 


A5)  eN-(r,t)  =  (QN/d)  +  (QN/2q)  S 


e-iu)t  + 


E  p2+a2 

^  W(e, -U), -T])  ) 

^  %  -K-a  ^ +  W(e. -T)))  j  + 

■  ^e^K  +  W(e,0,-T1)) 


where 


A6)  W(E,a),Tl)  =  g 


■F  ®k+£"®k  "^('"+i'n) 

We  again  view  the  potential  energy  of  the  ion  at  the  origin  as  heihg 
made  up  of  two  parts :  a  direct  interaction  potential  from  all  the  other  ■ 
ions,  £ind  a  potential  energy  of  interaction  with  the  displaced  electron 
charge  dlstrihution.  The  potential  energy  of  this  ion  is  given  hy 

A7)  $(rj  =  J*  QVj(r-rij^)  YSjj^QNCr-^-Rjjj)  +  eNg(r,t)  )  d^r 

n  '  ' 

where  QVj(r)  and  QN(r)  are,  respectively,  the  electrostatic  potential 
and  charge  density  of  an  ion  sitmted  at  the  origin.  For  point  ions, 
we  find  for  the  force  acting  on  the  ion  at  the  origin 
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p2+8i,2 


A8)  (W^N/n)  2  p  —  £  (£.e)  "  ^£,K+a^  ■" 

[  ^£,  K  ^  ^  ^ )  cosujt  + 

Re  [¥(£,  00,7])  +  W(£, -00,-1])  }  cbsoot  + 

.1 


-  6 


E>k+3 

■  ^P,K+£  [W(e, 00,-1])  +  ¥(£,-00,-7])}  Sinu)tj 

The  force  is  seen  to  be  somewhat  out  of  phase  with  The  sinoot  term 


in  Fj^  will  cause  ener^  dissipation,  and  implies  that  (Al)  does  not 
solve  the  equations  of  motion  exactly  for  the  lattice.  But  if  the  energy 
dissipation  is  slow  enough,  conpared  to  oo,  (Al)  will  still  give  a  nearly 
exact  solution.  In  this  case,  which  is  the  case  of  interest,  we  can 
identify  the  cosoot  term  in  F^  driving  force  on  the  ion,  and  the 

sinoot  term  as  a  dissipative  force.  Frequencies  are  then  calculated  as 
in  Part  I.  The  rate  of  energy  dissipation  per  unit  energy  is 
^  Ptt/oo 

A9)  1/t  =  (l/Mjo^)(oo/2n)  J  ~ 

o 

o  2 

=  (J4-TrqrN/2Maon)  ^£  2(  2+  2)  Im [w(£,  00,7] )+¥(£, -u), -7])  } 

The  expressions  for  frequency  and  decay  rate  can  be  combined  into 
the  expression 
AlO)  00^  +  2ioo/T  =  B(£,e,(o) 
where 

All)  B(^e,oo)  =  (W^Pn/MTj)  S  -  (£'e)^  X 

4TTe2 


X  lim  ^ 
7]  -*  0 


-  6. 


'"p,K+qT^  ■  ~^~~2 

K  [i  -  ^5^  (w(£,0,7])+¥(£,0,-7)))j 
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1/t  is  to  te  Identified  as  the  phonon  decay  rate.  The  equations  (AIO) 
and  (All)  determine  (u  im^plicitly,  but  we  will  find  that,  for  frequencies 
of  ordinary  magnitude,  the  real  part  of  W(£,u),Tl)  is  to  a  very  good 
approximation  independent  of  cu. 

We  shoijld  note  here  that  the  derivation  of  Eq.  (AIO)  has  been  more 
specialized  than  is  necessary.  In  general,  the  time  dependent  electron 
density  in  the  vibrating  lattice  can  be  written 
A12)  N(r,t)  -  S  K  -  + 

+  n(K+^(B)ei  [(K+£).r-a)t]  +  /(K+^c„)e-i  [(K+^)-r-u)t]^ 

This  form  follows  from  the  reality  of  N(r,t)  and  from  the  fact  that  the 
electron  gas  can  only  take  up  crystal  momentom  from  the  vibrating 
lattice.  Using  this  form  for  the  electron  density,  we  arrive  at  a 
relation  of  the  form  of  Eq.  (AlO). 

b.  Properties  of  the  Function  W(jo,io,T|) 

In  the  example  treated  in  Sec.  a,  the  real  and  imaginary  parts  of 
the  function  B(cL,e,u))  can  be  found  by  investigating  the  analytic  struc¬ 
ture  of  the  function 

•AI3)  W(£,cu'*')  =  lim  W(£,u),T|) 

'll  --  0 

where  the  notation  ou"^  means  that  u)  is  taken  to  have  a  vanishingly  small 
positive  imaginary  paid;.  W(p,(B'‘')  is  defined  by  Eq.  (a6)  as  a  sum  over 
energy  denominators  j  converting  the  smnmatlon  in  (a6)  into  an  integration 
gives 

Al4)  w(£,(i)+)  =  (STT^n)"^  J*  ^(fi/m*)(k*£  +  - 
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I 

I 


where  the  Integration  runs  over  all  k  suc5h  that  k<kj,.  This  integral, 
considered  as  a-  function,  of  coii5)lex  o)  (parametrized  hy  £),  has.  a  "branch 
cut  on  the  real  axis  of  the  complex  oj-plane.  The  imaginary  part  of 
W(£,u)+)  is  just 

AI5)  Im  [w(£,ou''')}  =  (8tt^)  ^  J  6  [(h/m)(k.£  +  -co]  d^k 

This  integration  can  be  performed  simply,  but  in  order  to  clarify 
the  general  featvires  of  the  problem  we  picture  the  integration  in  Pig. 

Ai.  The  dashed  line  in  Fig.  Al  represents  the  surface  K  in  k-space 
where,  for  fixed  £  and  co,  the  energy  denominator  in  the  integrand  of 
Eq.  (a14)  vanishes.  The  integration  for  the  imaginary  part  of  W(£,u)’*') 
just  runs  over  the  area  intersected  by  the  Fermi  surface  on  the  surface 
R,  and  is  proportional  to  this  area.  Now  since  the  energy  denominator 
is  a  linear  function  of  k  and  ou,  the  position  of  the  surface  R  travels 
linearly  with  ou  as  cu  is  varied.  We  call  u)©  a  value  of  co  for  which  R 
becomes  tangent  to  the  Fermi  surface.  Then  as  ou  -♦  coq  in  such  a  way 
that  R  goes  from  inside  the  Fermi  surface  to  become  tangent  to  it,"  the 
area  intersected  on  R  goes  to  zero  linearly  with  ou  -  ouq.  Thus  the 
imaginary  part  of  W(£,ou+)  goes  to  zero  linearly  as  ,ou  -*  coo  •  For  .values 
of  ou  for  ^ich  R  does  not  intersect  the  Fermi  surface,  ■w(£,ou’‘')  will  have 
no  imaginary  part. 

Performing  the  integration  for  Im  (w(£,u)+)}  gives 
A16)  Im[¥(£,ou+)}  =  -(m*3/8nh^p3)  [00  -  (ti/“*)p(ip-%) ]  [«>  "  (V“*)p(|p+%) ] 

...  for  (li/m*)p(|p-kp)<5o<  (li/in*)p(4p+l^) 

0  for  0U+  outside  this  range. 


i 

1 


I 


.) 
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which,  indeed,  falls  linearly  to  zero  as  co  approaches  the  values  for 
which  R  is  tangent  to  the  Fermi  surface. 

The  existence  of  am  imaginary  part  of  w(£,u)"'')  just  reflects  the 
existence  of  the  branch  cut  on  the  real  axis  of  the  cu-pleme.  From 
Eq.  (A14)  it  follows  that  W(£,u))  has  a  discontinuity  in  its  imaginary 
part  as.  ou  goes  from  above  to  below  the  branch  cut;  this  discontiniiity 
corresponds  to  a  change  in  sign  of  the  imaginary  part.  IHien  Im[w(£,a)+)} 

. can  be  seen  to  have  Just  half  the  value  of  this  discontinuity. 

From  inspection  of  Eq.  (A14)  we  see  that  W(£,u))  falls  off  as  l/ou 
for  large  ou,  and  that  the  only  non-ahalyticity  in  this  function  is  the 
blanch  cut  on  the  real  axis.  The  real  part  of  w(£,ou+)  can  then  be  found 
by  a  simple  contour  integration.  ¥e  can  write 

AI7)  f  — -•  ds  =  0  =  ds  -  inW(ji,u)+) 

where  C  is  the  contour  drawn  in  Pig.  A2.  Taking  the  imaginary  part  of 
Eq.  (AI7)  gives 

A18)  Re  [W(£,cu+)}  =  (1/tt)  Im[w(£,u)+)3 

s  -  0U+ 

The  integration  in  fact  runs  only  along  the  branch  cut  from  left  to 
right.  Performing, the  integration,  we  find 

AI9)  Re{;W(£,ou+)3  =  ((V“*)p%C(*/“*)p^- “  2ou]  + 

+  [ou  -  (*/m*)p(|p-kj,)][io  -  (h/m*)p(|p+kp) ]ln I 

I  (*/ m*)p(^+kji)  -0) 

It  should  be  noted  that,  for  small  u),  this  Just  becomes  the  function  W(£), 
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where  W(£)  is  the  function  which  has  been  found  to  he  responsible  for 
"kinks"  in  phonon  and  spin-wave  dispersion  curves.  The  nature  of  the 
anomolous  behavior  of  Re[w(£,u)+)}  in  the  regions  ou  ~  (ft/in*)p(ip-kp)  aiwi 
(ju  ~  (ili/m*)p(^+kj.)  is  determined  by  the  fact  that  Im{w(£,u)  )}  goes 
linearly  to  zero  in  these  regions.  For  u)  ~  (h/m*)p(ip-kji'),  for  instance, 
the  important  contribution  to  the  integral  (A18).  comes  from  the  region 
s  ~  (fi/m*)p(^-kp),  where  the  integrand  is  approximately 


A20) 


(h/m*)p(-^-kp) 

S  -  UU+ 


) 


0 


for  (ti/m*)p(aP-l^F)<s 
for  s  <(fi/“*)p(aP-%) 


and  this  region  contributes  a  term 

A21)  -K[uj  -  (h/“‘*)p(lP-kF)]ln|a)  -  (Vm*)p(ip-kF)  | 

to  Re  [W(£, («■'■)}. 

We  are  particularly  interested,  not  in  the  function  W(£>cu+),  but  in 
the  s\im  w(£,uH-)  +  W(jD, -u)+).  To  find  w(£, -(«+)  we  write 


A22)  r  ^ )  ds  =  0  =  -P  ^  )  ds  +  iTTW(£, -ou"^) 

C  s  +  0U+  s  +  (U+ 

where  C  is  the  contour  drawn  in  Fig.  A3,  and  s"  is  taken  to  have  a 
vanishingly  small  negative  imaginary  part.  Taking  real  and  imaginary 
parts  of  this  expression  yields 


A23)  lm{w(£,-a)+)}  =  (1/tt)  ds 

’  ■  "  s  +  0U+ 

Re  [W(£, -0)“^) }  =  (l/rr)  ^  ^[W(£,  s'^)} 

""  s  + 

where  advantage  has  been  taken  of  the  relation;  Im[w(£,s+)}  =  -Im{w(£,s~)}. 
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The  analogous  equations  for  W(£,u)+)  are,  t^lng  real  and  imaginary  parts 
of  Eq,.  (AIT), 

A2k)  Im{w(£,(u  )}  =  "(l/tr)  ds 

.  S  -  GU+ 

Re[w(£,c«  )}  =  (1/rr)  Iin[w(£,s+)}  ^3 

+■ 

S  -  U)+ 

Coinparison  of  (A24)  with  (A23)  yields  the  following  rule  for  constaniction 
of-  w(£,-u)+)  from  w(£,u)+):  Formally  replace  to  in  the  latter  function  by 
-to,  and  change  the  sign  of  the  imaginary  part.  It  follows  that 
Re[w(£,to+)  +  W(£, -io+)]  is  an  even  function  of  to,  whereas  Im[w(p,to+)+W(£, -u)+)} 
is  odd. 

We  wish  to  know  the  behavior  of  these  functions  as  functions  of  £ 
for  physical  values  of  the  frequency,  which  are,  in  order  of  magnitude, 
10"3(iik^/m*).  Then,  except  for  values  of  p  such  that  the  end  of  the  branch 
cut  falls  very  near  the  origin  in  the  co-plane,  we  can  expand  the  functions 
for  small  co  to  get 

.A25)  Im[w(£,cu+)  +  W(£, -u)+)}  =2(0:^  Im[V7(£,u)+)]j 

doo  CO  =  0 

Re[W(£,co+)  +  W(£, -CU+)}  =  2  Re[W(£,0)3 

These  relations  hold  for  all  values  of  p  except  those  for  which  the 

quantity  (fikp,/m*)(p-2kj,)  is  of  order  co.  But  since,  for  physlceil  values 

of  CO,  these  values  of  p  are  within  10” ^k^,  of  2kp„  measured  spectral 

properties  will  in  effect  be  given  by 


A26)  5lm{w(£,(o'*‘)  +  W(£, -(«■*■)} 


CO 


m*^/  BrfhSp 

0. 


for  p<2kp 
for  p>2kp 


iHe[W(£,co+)  +  W(£,-co+)}  =  W(£) 
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where  w(£)  is  the  function  already  fovind  to  he  responsible  for  "Mnks" 
in  phonon  dispersion  curves. 

A  short  summaiy  of  the  results  of  this  section  is  in  order.  We 
have  found  the  function  W(£,u))  to  have  a  branch,  cut  on  the  real  axis  of 
the  complex  (u-plane.  The  positions  of  tfie  end  points  of  this  branch  cut 
depend  parametrically  on  £.  Near  the  ends  of  the  branch  cut  Im {¥(£,(«■*■)} 
falls  off  linearly  to  zero,  and  remains  zero  off  the  branch  cut. 

By  contour  integration  in  the  (o-plane,  it  has  been  established  that 
Re[W(£,u)+)}  has  anomalies  of  the  form  ((i)-u)Q)ln|(i)-a)Q |,  where  ooq  is  the 
end  of  a  branch  cut.  It  has  also  been  established  in  this  fashion  that 
the  functions  of  physical  interest,  in  cases  where  frequencies  are  small, 
are  Re{w(j£,0)  }  and  U3  Im[W(£,u)+)}  (  ^  ,  and  that  these  functions, 

considered  as  functions  of  £,  display  exactly  the  behavior  previously 
fotmd  for  anomalies  in  phonon  dispersion  curves  and  phonon  lifetime  curves 

In  particular,  we  have  found  that,  for  p  ~  2kj',  the  branch  cut  has 
one  end  near  the  origin  of  i^e  oo-plane.  Near  this  end  of  the  branch  cut, 
we  have 

A30)  Im{w(£,u)+)}  K[u)  -  (*/in*)p(ip-kp)] 

It  then  follows  that  anomalies  in  the  phonon  decay  rate  take  the  form  and 
magnitude 


A31)  Im[W(£,(«+)  +  W(£,(u+)3 


K  0) 


for  p<2kji 


G  for  p>2kj' 

/ 

and  by  contour  integration  it  can  be  established  that  anomalies  in  the 
phonon  dispersion  curves  have  form  and  magnittide 
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A32)  -K  V  (p-2kp)  ln|p-2kii.| 

where  v  =  ftp/2m*  is  the  velocity  of  an  electron  whose  wavenumber  is 
c.  Non-Spherical  Fermi  Surfaces 

In  cases  where  the  Fermi  surface  is  not  spherical,  the  anomalous 
behavior  of  dispersion  curves  and  lifetimes  can  still  be  expected  to  ai;ise 
from  the  behavior  of  the  function 

A33)  =  (8tt3)-^  J 

;daere  the  integral  runs  over  all  electron  energy  states  beneath- the 
Fermi  surface.  This  integral  may  be  quite  difficult  to  evaluate  for 
general  Fermi  surfaces.  We  show  in  this  section,  however,  that  in  the 
im^rtant  region  near  the  end  of  the  branch  cut  in  the  cu-plane,  for 
small  (ju,  this  function  has  the  behavior  already  observed  in  the  case  of 
a  spherical  Fermi  surface. 

The  imaginary  part  of  W*(e,u)‘'^)  is  given  by 

A34-)  Im[W*(£,uj'^)}  =  (8n%)"'l  J  6 

This  integral  runs  over  the  surface  R  in  k  space  defined,  for  fixed  £ 

and  u),  by 

«5) 

The  behavior  of  Im[W*(£,- («■*■)},  in  the  important  regions  near  the  ends  of 
the  branch  cut  of  W^(e,u)),  is  determined  by  the  shape  of  the  surface  R 
for  values  of  cu  which  make  R  nearly  tangent  to  the  Fermi  surface.  Inspec¬ 
tion  of  Fig.  A^l-  indicates  that,  for  £  in  such  a  direction  that  R  becomes 
tangent  to  a  convex  part  of  the  Fermi  surface,  the  imaglnaiy  part  of 
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V^(£,u)+)  will  indeed  fall  off  linearly  to  zero  as  ou  approaches  the 
branch  point.  The  behavior  found  in  the  spherical  case  should  then  be 
repeated  in  this  case^  for  £  in  proper  directions. 

We  proceed  to  find  an  expression  for  W*’(p,(D+)  when  u)  is  near  the 
end  of  the  branch  cut,  for  such  £  that  the  region  of  tangency  occurs 
where  the  Fermi  surface  is  convex  and  for  small  co.  It  will  be  con¬ 
venient  to  consider  the  properties  of  the  constant  energy  surfaces  in 
the  £  direction,  since  we  shall  find  that  R  becomes  tangent  to  the  Fermi, 
surface  in  this  direction.  In  the  vicinity  of  the  point  as  indi¬ 
cated  in  Fig.  A5,  we  can  write  in  the  form* 

A36)  \  a(k^+ip^)  +  bCk^+ip^)^  +  cCk^+ip^)^ 

where  k^  is  the  component  of  k  normal  to  the  surface  Ej^  =  Ej^  at  the 
point  -^£,  and  the  directions  of  k^  and  k^  have  been  chosen'  'to  eliminate 
cross  terms  in  k^k^  from  Eq.  (A34).  We  can  likewise  write. 

-  a(k^+|p3^)  +  bCk^+lpg)^  +  c(k^+|p^)^ 

.  so  the  defining  equation  for  the  surface  R  becomes 
A38)  -2a(k^+^^)  -  fiou  =  0 

The  branch  cut  then  has  its  end  for  small  w  when  £  is  in  the  vicinity  of 
2kQ,  where  kg  is  the  vector  on  the  Fermi  surface  in  the  £  direction. 

For  such  £,  the  imaginary  part  of  W*(£,u)+)  is  given  by 
A39)  ImfW^ (£,(«+) 3  =<  (A/8rr%)  J  6  [(2a/h)(k3_44pi)  +  ou] 
where  A  is  the  area  intersected  by  the  Fermi  surface  on  the  surface  R. 


^Inversion  symmetry  has  been  assumed  of  the  Fermi  surface. 
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The  equation  of  the  Fermi  surface  near  the  point is 
AkO)  Ejj,  -  +  a(k^+^^)  +  'bCk^+iP^)^  +  cCk^+ip^)^ 

The  area  intersected  on  R  is  thus  elliptical  in  sha;^,  ^d  the  area  A  is 

1 

A41)  a  =  tt(T3c)“^  (iL,  -  El  +  = 

a 

-  ^TTft(bc)  ^  [(2a  cos0)(|p-k^)/h  +  (ju] 

where  0  is  the  angle  between  the  k^^  and  £  directions.  Then  the  ima^nary 
part  of  W^(£,UD+)  is 

A42)  Im[W*(£,a)+)}  =  |  ggTfa^b^'l '  cos0)(|p-k^)/h] 

—  for  (2a  cos0)(^-k^)/h  <  U) 

0  for  u)  <  (2a  oos0)(-|p-k^)/ft 

The  constants  a,  b  and  c  have  simple  physical  interpretation.  The 
ratios  (a/2b)  and  (a/2c)  are,  respectively,  the  radii  of  curvature  and 
R^  of  the  surface  Ej^=E^  at  the  points  ±  •^.  The  constant  a  can  be 
Interpreted  by  noting  that 


A43)  V  =  \  \  1  =  S 

-  -  k=-|£ 

where  v  is  the  velocity  of  an  electron  whose  wave  vector  is  The 

vector  n  is  the  unit  vector  normal  to  the  surface  E,  =  Ei  at  the  point 
—  k 

Prom  Eq.  (a42),  using  the  results  of  Sec.  b,  we  see  that  the  phonon 
decay  rate  can  be  expected  to  display  anomalies  of  form  eind  magnitude 


A41^)  Im{;w*'(p,u)+)  +  W^(£,-u)+)3  = 


K  u) 


0 


where 
Ai^5)  K  = 


for  p<2k^ 
o 

for  p>2k^ 
o 


eind  the  anomalies  of  the  phonon  dispersion  curves  will  have  form  aM 
magnitude 

Ait6)  -K  V  COS0  (p-2k^)  ln|p-2k^|  . 
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APPENDIX  2 

The  Electron- Ion  Matrix  Element  in  Aluminum 

1 

Bardeen's  calcvilatiaa  of  the  matrix  element  (k+£|vv(r)  |k)  uses 

the  basic  assvimption  that  the  electron  wave  functions,  -written  in  the 

Bloch  form  e^— *—  ,  can  be  approximated  as  u^Cr)  e^— ,  where 

Uo(r)  is  the  lowest  eigenfunction  of  the  conduction  electron  Hamiltonian. 

This  assumption  is,  pres\niiably,  -valid  for  mono-valent  me-bals;  we  shall 

attempt  to  show  here  that  it  can  be  replaced  by  two  other  assunptions  ; 

that  the  electron  wa-ve  functions  can  be  adequately  represented  by  single 

IT 

Orthogonalized  Plane  Wa-ves  (OPW),  '  and  that  each  ion  core  is  well  con^ 
fined  -wl-thin  its  own  lanit  cell. 

This  replacement  of  Bardeen's  assumption  is  possible  because  -the 
assumption  has  been  used  only  to  determine  the  properties  of  in 

the  -vicinity  of  the  unit  cell  wall.  The  impor-bant  properties  6iie, 
first,  that  u^(£)  should  ha-ve  spherical  symmetry  in  this  region,  and 
second,  -that  shpuld  not  -vary  much  with  k.  The  OPW  are  constructed 

by  adding  to  plane  waves  linear  coiubinations  of  ion  core  wavefunctions 
(wi-th  coefficients  so  chosen  as  to  sake  the  OPW  orthogonal  to  all  core 
s-tates.)  For  ions  which  are  well  confined  within  their  own  unit  cells, 
■then,  the  OPW  necessarily  have  almost  pure  plane  wa-ve  character  near  the 
cell  wall.  Then  Uj^(£)  -will  have  spherical  symmetry  in  this  region,  and 
will  in  fact  be  very  nearly  flat.  The  -value  -will  depend  on  k 

only  -through  a  normalization  cons-tant.  If  the  largest  ion  core 
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wavefunction  falls  off  as  e  '  ,  then  the  natural  parameter  for  a  power 

?  l8 

series  expansion  of  the  normalization  constant  is  (ka)~.  From  Froese's 

calculations  for  the  trivalent  aluminum  ion,  for  instance,  we  find 

kj,a  =  1/  5,  so  the  variation  of  ujj.(rg)  will  he  small  for  k  in  the  region 

of  interest,  (it  is  6l1so  desirable,  though  of  less  importance,  that 

the  nxamber  C  in  Eq.  (l,23)  he  independent  of  k.  C  is  proportional  to 

2  2 

the  second  derivative  3  u^(r)/3r^,  evaluated  at  the  cell  wall.  This 

quantity  will  he  quite  small  if  the  single  OPW  assumption  is  justified, 

since  it  falls  off  with  distance  in  the  same  way  as  the  core  wavefunctions 

If  a  power  series  expansion  is  made,  the  parameter  is  ka  =  k/5kp,  so  the 

variation  with  k  should  he  reasonably  slow. ) 

Bardeen's  calculation  of  the  electron-ion  matrix  element  therefore 

goes  through  under  our  assumptions.  The  question  of.  whether  the  single- 

OPW  approximation  is  good  for  any  particialar  metal  is  rather  difficult 

17 

to  answer,  hut  the  work  of  Heine  shows  rapid  convergence,  in  'aluminum, 
for  the  expansion  of  the  Bloch  functions  in  single  OPW's. 
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APPENDIX  3 

Competing  Decay  Processes  for  Phonons 


a.  Three -Phonon  Effects 

Terms  in  the  potential  energy  of  the  lattice  which  are  cubic  in 

the  ionic  displacements  give  rise  to  phonon-phonon  interactions.  These 

interactions  contribute  to  the  phonon  decay  rate,  and  can  be  expected  to 

reduce,  to  some  extent,  the  relative  importance  of  electron-phonon 

effects.  We  show  here  that  this  reduction  is  not  serious,  estimating 

the  decay  rate  from  phonon-phonon  interactions  by  using  a  phenomeno- 

19 

logical  treatment  due  to  KLemens. 

The  cubic  terms  in  the  lattice  Hamiltonian  can  be  written  in  second 
quantization  as  terms  cubic  in  the  phonon  field  operators.  If  only  one 
phonon  is  initially  present  in  the  lattice,  having  wave  number  and 
polarization  e,  the  only  cubic  terns  in  the  Hamiltonian  capable  of 
destroying  the  phonon  are 


CteE-2;£.i.£) 


where  a  and  a  are,  respectively,  the  operators  which  create  and 

destroy  phonons  of  wavenumber  ^  and  polarization  e.  These  operators  have 
19 

the  matrix  elements 


A48)  (n'l  a  I  n)  =  [hn/2Mu]^  ^n',n-l 

(n'l  a*l  n)  =  [h(n+l)/2Hi)]2 

where  the  polarization  and  wavenumber  subscripts  have  been  suppressed. 


I 


4o 


We  follow  Klemens  in  taking 

A49)  C(k,£,£;a,t,e  )  =  S  ^ 


(3N3)t 


where  y  ~  2,  v  is  the  phase  velocity  of  the  phonon  initially  present, 
and  the  sum  runs  over  the  Bravais  lattice. 

The  transition  probability  for  decay  of  one  phonon  into  all  possible 
pairs  of  phonons  is,  in  the  Debye  approximation, 

A50)  W  =  E3  .  (c^c^/Cg)  q  J  k|k-K-a|  eChc^k  +  hc^  |k-X-q( -ftc^q)  d^k 

where  the  integral  goes  over  the  Brillouin  zone,  and  Cg  is  the  sound 
velocity  of  the  phonon  whose  polarization  is  s.  The  integral  in  Eq.  (A50) 
can  be  done  rather  conveniently,  if  the  Brillouin  zone  is  taken  to  be  a 
sphere  of  radius  q^,  in  the  coordinates  u,  v,  cp,  where^*^ 

A51)  u  =  i(k  +|k-K-sD 

V  =  i(k  -|k-K-_aD 

and  cp  is  the  azimuthal  angle  about  K+£.  In  these  coordinates  the  transi¬ 
tion  probability  becomes 
Y^tm 

"  \,t 

where 

A53)  a  =  Cg  +  c^  ;  b  =  Cg  -  c^ 

and  the  integral  runs  over  the  area  of  the  u,v  plane  which  is  shaded  in 
Fig.  a6.  The  6-function  insures  that  in  fact  this  integral  shall  go  over 
only  the  line  drawn  dashed  in  Fig.  a6. 


A52)  W  =  y  ^  q  S  (c  c,/c  )  E  ^  ff  (u^-v^)^  6(au  +  bv  -cq)  du  dv 
'  ftirMW  ^  =  +•.  '  s  V  e '  K  Trr. — n-'  J  '  '  ' 


1 


kl 

The  two  polarization  "branches  require  six  different  decay  processes 
to  be  distinguished,  which  we  investigate  individually-  Each  of  these 

I 

processes  is  characterized  by  values  of  c/a  and  b/a.  Since  the  function 
■^|K+^|  has  the  value  -gq  for  K=0  and  minimum  value  gK-gq  for  ^^0,  it  is 
possible  to  know  from  Fig.  A6  and  the  values  of  c/a  and  b/a,  whether  for 
given  K  a  process  makes  any  contribution  to  the  transition  probability. 

We ' take  the  longitudinal  sound  velocity  to  be  twice  the  transverse  sound 
velocity,  and  we  symbolize  by  1  ->  H-t  the  process  in  which  a  longitudinal 
phonon  decays  into  one  longitudinal  and  one  transverse  phonon.  The  six 
decay  processes  are 

i)  t  -«  1+1  Here  the  dashed  line  in  Fig.  A6  lies  entirely  outside  the 
area  of  integration,  since  cq/a=-^  and  b/a=0.  This  result  is  independent 
of  the  value  of  K,  so  this  process  makes  no  contribution. 

ii)  t  -*  1+t  Here  again  the  dashed  line  falls  outside  the  area  of  inte¬ 
gration  for  all  K,  since  cq/a=q/3  .and  b/a=l/3* 

iii)  t  -»  t+t  Here  the  dashed  line  is  the  line  u=^.  For  K=0  we  get  a 
contribution,  but  for  k/O  there  is  none. 

iv)  1  -»  1+1  The  considerations  here  are  the  same  as  in  process  (ill). 

v)  1  -♦  1+t  Here  the  dashed  line  falls  inside  the  area  of  integration 
for  K=0,  and,  with  soire  values  of  q,  also  for  K  one  of  the  shortest 
reciprocal  lattice  vectors  (of  magnitude  2qo.)  The  values  of  q  which 
give  contributions  for  k/o  can  be  seen  to  lie  outside  the  sphere  of  radius 
2qo/3^  since  cq/a=T2q/3  and  b/a=l/3. 

vl)  1  -»  t+t  This  process  also  contributes  for  K=0  and,  with  q>2qQ/3^ 
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for  K  the  shortest  reciprocal  lattice  vectors,  since  cq/a=q.  and  'b/a.=0. 

In  view  of  the  above  survey,  we  coD^pute  the  transition  probability 
for  q<2qQ/3^  where  only  processes  (iil),  (iv),  (v),  and  (vi)  contribute, 
and  these  only  for  K=0.  Except  in  process  (v)  both  final  phonons  lie  on 
the  same  polarization  branch,  and  the  integral  in  Eq.  (A52)  becomes 

A54)  J'J(u^-v^)^6(u-cq/a)  du  dv  =  q5  [(c/a)^  -(c/a)^/6  +  l/80] 

for  ^<20^/3.  For  processes  (iii)  and  (iv)  c/a=^,  and  the  integral  takes 
the  value  qVsO*  For  process  (vi)  c/a  =  1,  and  the  integral  takes  on  a 
value  twenty-five  times  as  large.  In  process  (iv),  the  only  remaining 
contributor,  the  integral  in  Eq.  (A52)  becomes 

A55)  fj’(u^-v^)^6(u+bv/a-cq/a)  du  dv  ~  qV^O  • 

The  most  important  contribution  apparently  comes  from  process  (vi)  and, 
including  only  this  process  we  get,  for  constants  corresponding  to 
aliuninum, 

^  I 

A56)  W/ci)  0.06  (q/qo)^ 

Reference  to  Fig.  4  shows  that  this  decay  rate  is  quite  small  compared 
to  the  decay  rate  from  electron-phonon  interaction,  becoming  conparable 
only  near  the  boundary  of  the  Brillouin  zone. 

For  transverse  phonons,  the  only  decay  mechanism  is  process  (iil)# 
so  the  decay  rate  will  be  twenty- five  times  smaller.  ^ 

b.  Two-Pair  Creation 

One  possibly  important  phonon  decay  process  can  be  described  as 
follows:  the  phonon  excites  an  electron-hole  pair,  and  the  excited 
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electron,  without  falling  back  beneath  the  Fermi  sea,  subsequently 

excites  a  second  pair  through  electron-electron  Interaction.  Since 

energy  conservation  is  not  necessary  in  the  intermediate  state,  the 

restrictions  which  might  cause  a  break  for  phonons  of  wave  number  2%  do 

not  exist;  so  if  this  process  contributes  significantly  to  the  phonon 

decay  rate  it  will  partially  mask  out  the  effects  which  interest  us. 

We-  estimate  the  decay  rate  for  this  process  here,  however,  and  show 

that  it  is  much  smaller  than  the  decay  rate  into  single  electron-hole 
* 

pairs . 

We  reconsider  shortly,  first,  the  decay  rate  from  single  pair  crea¬ 
tion.  Ihis  decay  rate  can  be  calcxilated  as  a  transition  probability  for 
the  whole  system,  from  the  state  with  the  electron  gas  in  its  ground 
state  and  one  phonon  present  to  the  state  with  one  electron-hole  pair 
excited  and  the  phonon  absent.  We  represent  this  process  pictorially  in 
Fig.  A7,  taking  account  of  momentum  conservation.  Ihe  transition  proba¬ 
bility  is  given  by 


A57)  \  =  (l/8TT3h)(N/n)E  Jd^k  |Vp(^+K)|2  \ 

where  y^(£+K)  is  the  matrix  element  of  the  electron-phonon  interaction 
between  the  initial  and  final  states,  and  is  given  by 

A58)  V  (_a)  =  (e.£)  (h/2M«)2 


*Such  a  two-step  decay  can  occur  by  several  different  processes,  each 
of  which  will  contribute  to  the  two-pair  decay  rate  a  term  of  the  same 
order  as  the  term  estimated  here. 


1/  ’ 


where  the  Thomas -Fermi  approximation  has  heen  used  for  the  effect  of 
electron  shielding.  The  integration  on  k  is  restricted,  since  the  final 
state  must  have  an  electron  and  a  hole,  to  the  region  of  k  space  defined 

A59)  k<kj. 

|k+£+Kj>kp 

and  since  the  energy  conserving  6-function  can  he  written 

A60)  \  =  (m*/ft2)§[-k.(_^+K)  +  i(^+K)2  - 

the  region  available  for  the  integration  on  k  is  just  that  indicated  in 

Fig.  A8,  which  should  he  imagined  as  rotated  about  the  vector  £+K.  The 

unshaded  volixme  is  forbidden  by  the  electron-hole  requirement,  and  energy 

conservation  requires  that  k  lie  in  a  plane  perpendiciilar  to  ^+K,  and  a 

distance  -i|£+K|  +  m*(u/h(c^+K|  along  the  2+K  direction.  (Since  m*co/Kk^lO 

this  distance  is  almost  -■i|jl+K[.)  The  plane  is  indicated  by  the  dashed 

line  in  Fig.  A8.  k  then  lies  on  a  ring,  whose  area  is 

A6i)  a  =  n[k|  -  (i|^+K|  -  mVli|a+K|)^]-iT[k|  -  (i|£+K|  +  mV*|ja+K|)^]  = 

=  2TTm*u)/h 

The  Integration  over  the  6-function  introduces  a  factor  |£+K|~^.  The 
result  for  W^,  then,  estimating  that  about  10  reciprocal  lattice  vectors 
enter  into  the  sum  and  giving  |£+K|  the  value  2kp,  is 

A62)  =“  (5NArr^n)  |Vp(2kp)  f  =-  lO'^cu 

The  two -pair  creation  is  Illustrated  in  Fig.  A9,  with  the  electron- 
electron  interaction  represented  by  a  dashed  line  carrying  momentum  lin. 


r  . 

'4 


} 


^5 

The  transition  probability  for  this  process  can  be  written 

JS63)  w,  =  (V8n3)  jAfAfdSj,  llm  |Vj,(fl+K)P|v,(s)|2  ^ 

'll-  0  - - 

*  Vs+K  -  Ek-li(urHTl)  | 


where  the  pole  of  the  Integrand  has  been  treated  as  directed  by  Schiff 
^e^s)  matrix  element  of  the  (screened)  electron-electron  inter¬ 
action  between  intermediate  and  final  states,  given  by 

A6k)  vjn)  . 
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n^+ap- 


;  VQ(2kp)  =“  3  Vp(2kp)  for  aluminum 


for  Thomas -Fermi  screening. 

Requiring  that  the  final  state  contain  two  electron-hole  pairs 
restricts  the  sums  on  k  and  £  'by 


A65)  k<kp 

|k+£+K-n|  >kj, 


P<kp 
|£+nl  >kp, 


The  energy  differences  E,  ,  -E,  and  E  -E  must  both  be  positive,  so 

K+£+k-n  k  £+n  £  > 

each  must  be  smaller  than  or  equal  to  liu).  We  take  advantage  of  this 
fact  by  writing 

a66)  "  ■ 


=  5[k*  (£+K-n)+|-(£+K-n)^  -m^s/iljeEE'n+gn^  -.m**^(u)-s)/li]  ds 

where  the  limits  on  the  s  integration  are  required  by  the  necessity  for 
positive  energy  differences. 


The  regions  available  for  the  £  and  k  integrations  are  shown  in 
Fig.  AlO.  The  area  allowed  for  £,  for  fixed  s,  is  just  2Tmi*((u-s)/h^ 
and  a  factor  n  ^  is  introduced  by  integration  over  the  6-fuhctibn.  Per 
forming  the  £  and  s  integrations  then  gives  a  factor 


a67) 


(2u/n)  (_g+K-n)  +  i(£+K-n)^]}  for  0<(li/m* ) [k*  ( C|+K-n)+5(5[+K-n) 


|1)  otherwise 

The  k  integration  is  thus  confined  to  the  ring,  of  triangular  cross  sec¬ 
tion,  indicated  by  the  blackened  area  in  Fig.  AlO. 

The  reason  for  the  smallness  of  has  now  become  apparent.  Barring 
bad  behavior  of  the  squared  energy  denominator  in  Eq.  (A63),  integration 

O 

over  k  wotild  give  a  momentum- space  factor  u)'^,  and  since  only  phonon 

energies  and  electron  energies  enter  the  problem,  we  would  expect  to 

be  of  the  order  of  magnitude  (■hu)/E„)'^^a).  In  fact  we  now  investigate  the 

J? 

behavior  of  the  energy  denominator,  and  show  that  it  can  be  expected  to 

,  BO  that  Wg  still  remains, 

small  in  comparison  to  W^. 

We  can  write  the  squared  energy  denominator  in  Eq.  (a63)  as 

A68)  _ 1 _  _  m*  / _ 1 _ 

1^  2iT|n3  (^+K)+i(£+K)^-m  ((«+lTl)/li 

_ 1 _ 

k*  (£+K)+5(_a+K)^-m* (ou-iTl)/K 

This  denominator  behaves  badly  only  for  values  of  n  for  which  the  allowed 

region  of  the  k  integration  intersects  the  plane.  Indicated  by  a  dotted 

line  in  Fig.  AlO,  on  which  the  energy  denominator  becomes  zero.  For 


introduce,  at  worst,  a  factor  of  order  (Ej/hu)) 


hi 


such  values  of  n,  neglecting  for  the  moment  the  shape  of  the  integration 
volume  and  the  variation  of  factors  in  the  numerator,  integration  over 
the  coii5)onent  of  k  parallel  to  £+K,  gives 

)  IS,+K|  -  ni*(u)+iTl)/ft  \  ^ 


a69)  (m*/2iT^3|^+K|) 


In 


In 


m*(uH-iTl)/Bi  ' 

V  m*(uH-iTl)/l5  / 


since  the  minimum  value  of  k^^  is  -i'|£+Kj .  There  is  no  contribution  from 
either  limit  from  the  real  parts  of  the  logarithms.  The  imagimry  peirts 
of  the  logarithms  contribute,  in  the  limit  T|  -►  0,  a  factor 


A70)  (m*/h3|q+K|)  [l/uj  -m*/  C^la+K| (k.,^^+ila+K| ) ]} 

the  second  term  being  negligible  in  comparison  with  the  first.  The 
simplest  procedure  available  to  us  at  this  point,  then,  will  be  to  over¬ 
estimate  the  effect  of  the  squared  energy  denominator  by  giving  it  the 
value  (m*/b3|£+K|2(ju). 

The  allowed  volume  for  the  k  integration  is  just  m*^uj^/li^  |q+K-n|  > 
and  the  function  given  in  (a67)  has  the  average  value  n  co/n  over  this 
volume.  The  integration  on  n  can  be  approximated  by  giving  n  and  j^+K-nJ 
the  value  2kp,  and  since  the  n  integration  cuts  off  at  n=2kp  the  integra¬ 
tion  volxime  is  32TTkp/3.  The  range  of  the  reciprocal  lattice  sum  here  is 
|£+K|<4kp,  so  about  80  -terms  can  be  expected  to  enter  into  -this  sum. 

The  final  result  is 

m*3(u  .o  _l, 

A71)  |ve(2kj.)  r  10 


/ 


so  the  decay  rate  from  two -pair  creation  can  be  neglected  in  coni«riBOn 
to  the  decay  rate  from  single  pair  creation. 
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i^PEKDIX  k 

Images  in  Spin -Wave  Spectra  of  Antiferromagnets 

The  negative  spinrwave  "dispersion  curve"  of  Fig.  6  implies  that> 
for  the  model  used,  the  assvuiption  of  a  ferromagnetic  ground:  state  is 
not  justified.  Rather,  it  implies  that  the  ground  state  of  the  spin 
system  sho\ild  he  an  antiferromagnetic  spiral  structure,  whose  pitch  c^ 
he  deduced  from  the  position  of  the  minimum  in  the  negative  dispersion 
curve.  ^  The  extreme  sensitivity  of  our  results  to  dmall  changes 

in  the  position  of  the  Fermi  siirface  asstires  that  these  comments  should 
not  he  assigned  more  than  a  speculative  significance,  hut  it  is  inter¬ 
esting  to  note  that  the  rare  earths  do  display  spiral  structure  at  some?- 

24 

what  elevated  tenperatures.  It  may  he  that  the  sensitivity  of  the 

dispersion  curves  (suggested,  in  fact,  hy  the  disparity  between  the  two 

c\irves  of  Fig.  6)  can  help  to  explain  some  of  the  reinarkahle  magnetic 

24 

properties  of  the  rare  earths. 

The  results  of  Part  II  are  applicable  to  rare  earth  metals  only 

when  in  a  ferromagnetic  state,  either  naturally,  or  under  the  Influence 

24 

of  a  strong  magnetic  field.  For  spiral  antiferromagnetic  spin  arreuige- 
ments  on  Bravais  lattices,  Kasuya?’^  as  well  as  Yosida^^  has  found  a  spin 
wave  spectrum  given  hy 

AT2)  [Cd  - 1(^)  +  I(^)][i(s)  -  i  Kso+s)  -  5  Kao-a)]}^ 

where  d  is  a  constant  used  to  describe  anisotropy  energy  in  the  hard 
direction  of  nagnetizatlon,  and  l(^)  is  proportional  to  our  2jj(w(^K)-W(K)). 


50 


;3o  gives  the  pitch  of  the  ground-state  spiral  structure.  Eq.  (A59) 

indicates  that  in  the  spin-wave  spectra  of  antiferromagnetic  stjructvires 

images  of  the  Fermi  surface  will  appear,  with  magnitude  comparahle  to  | 

the  magnitude  they  have  for  ferromagnets.  The  condition  for  a  kink  is  I 

I 

that  £  ^tisfy  one  of  the  conditions  j 

A73)  1^  K|  =  2kp 

h  t  ao+  K|  =  2ky  . 
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FIGURE  2 

Plot  of  the  function  w(£).  Note  infinite  slope  at  p=2kj,  . 
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FIGURE  3a  (previous  page) 

Dispersion  curves  calculated  from  the  Uniform  Electron  Gas  Model  for  con¬ 
stants  corresponding  to  aluminiun.  Ihese  curves  are  for  the  (ill)  direc¬ 
tion  of  reciprocal  space  j  "l"  and  "t"  refer  to  longitudinal  and  trans¬ 
verse  directions  of  polarization.  Fermi  surface  images  are  indicated  ■by- 
arrows.  While  the  considerations  of  Sec.  B  suggest  that  the  model  used 
is  inappropria-fce  for  aluminum,  kinks  of  this  general, order  may  well 
appear  in  other  cases. 

FIGURE  3t>  (a'bove) 

Characteristic  "upward"  and  "downward"  anomalies  of  the  dispersion  CTirves. 
Scales  are  enlarged  five  times  from  the  scales  of  Fig.  3a*  While  no 
simple  downward  anomaly  occurs  in  the  case  of  aluminum,  the  example  shown 
here  should  be  helpful  in  in-fcerpre-tation  of  such  einomalles  where  they  do 
appear,  as  in  the  case  of  lead.  (Fig.  If) 


FIGURE  4 

Measured  longitudinal  'brauKdi  of  the  phonon  dispersion  curve  for  ‘the 
(no)  direction  in  lead.  (After  Brockhouse  ^  al;  see  Ref.  4.)  The 
qualitative  similarity  of  this  curve  to  the  theoretical  curves  of 
Fig.  3  is  striking.  Cie  anomaly  rcughly  indicated  by  the  arrow  is  a 
"downward"  type. 
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FIGURE  5 

Calculated  phonon  decay  rate,  expressed  as  a  fraction  of  frequency,  for 
£  in  the  (ill)  direction,  "x"  and  "t"  indicate  longitudinal  and  trans¬ 
verse  polarization  tranches.  The  Uniform  Electron  Gas  model  has  heen 
used  with  constants  corresponding  to  aluminum.  Decay  rates  for  small 
£  are  controlled  hy  more,  conplicated  processes,  (see  Ref.  lO)  and  are 
'  not  given  correctly  hy  these  curves.  Again,  the  model  used  is  probably 
inappropriate  for  aluminum,  but  such  breedcs  are  expected  to  appear  in 


other  cases. 


WAVENUMBER 


2  k, 


FIGURE  6 

Calculated  spin-wave  dispersion  curves  In  hig^-symmetry  directions. 
Fermi  surface*  images  are  indicated  ty  arrows.  Frequency  values  are 
calculated  using  the  coupling  constant  G  for  gadolinium,  derived  from 
the  observed  resistivity;  m*  is  taken  to  be  the  free-electron  mss. 
Frequencies  for  other  rare-earth  metals  can  be  obtained  by  multiplying 
by  (g-l)j/7.5  . 


Region  of  the  k  integration  for  W(£,(i)+) 
for  fixed  £  az3d  cu,  the  aurface  i^re  th 
(iU.4)  vanishes. 


The  dashed  line  represents 
energy  denominator  of  Eg,. 


PIC5URE  A5 


Approximation  of  the  properties  of  the  constant  energy  svirfaces  for 
evaluation  of  W* (£,(«■*■),  for  a)  near  the  end  of  the  branch  cut  and  £ 
approximately  a  "diameter"  of  the  Fernd.  surface.  Hhe  surface  R  is 
seen  to  lie  in  the  region  k  =“  . 
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